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We investigate the interaction of trains of femtosecond microjoule laser pulses with dielectric materials by 
means of a multi-scale model. Our theoretical predictions are directly confronted with experimental obser¬ 
vations in soda-lime glass. We show that due to the low heat conductivity, a significant fraction of the laser 
energy can be accumulated in the absorption region. Depending on the pulse repetition rate, the material 
can be heated to high temperatures even though the single pulse energy is too low to induce a significant 
material modification. Regions heated above the glass transition temperature in our simulations correspond 
very well to zones of permanent material modifications observed in the experiments. 


The use of femtosecond (fs) laser pulses for modifica¬ 
tions of transparent materials is nowadays a well estab¬ 
lished technique with many applications in science and 
engineering. Material modification caused by a pulsed 
fs laser, i.e. a train of fs pulses, can be controlled by 
numerous experimental parameters such as repetition 
rate (RR)^, incident energy^, or focusing conditions^. 
Thus, ultrashort intense laser pulses are a versatile, 
highly adaptable processing tool for micro-machining of 
various materials^, including human tissues in the con¬ 
text of ophthalmic surgery^. 

When a fs laser pulse is focused into a dielectric ma¬ 
terial to intensities exceeding a few TW/cm^, a consid¬ 
erable number of electrons undergo photo- and impact- 
ionization processes^: electrons from the valence band 
(VB) are promoted to the conduction band (CB)^. Once 
the carrier density in the CB becomes non-negligible, 
these electrons strongly influence the pulse propagation 
dynamics, and the laser energy is efficiently absorbed in 
the near focal-region^. At the end of the irradiation pro¬ 
cess, typically after a hundred fs, the electrons in the CB 
transfer their energy to the lattice through collisional pro¬ 
cesses, leading to an increase of the material temperature 
in the focal volume^. On longer time scales, this energy is 
transferred toward the surrounding cold matter through 
heat diffusion^^. Provided that the amount of laser en¬ 
ergy absorbed in the focal volume is sufficiently large, 
permanent material modifications are induced. This can 
be achieved either with a single fs micro joule laser pulse, 
or by using cumulative heating techniques'^, i.e., multi¬ 
ple pulses or pulse trains. Due to the low heat diffusivity 
in dielectric materials 10“^ cm^/s), the cooling time 
over the typical size of the focal spot (a few microns) is 
of the order of ^ 10 ps. Therefore, by adjusting the time 
interval between consecutive laser pulses appropriately 
(kHz-MHz RR), one may accumulate the laser energy 
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in the absorption region and reach very high tempera¬ 
tures, even if the single pulse energy is too low to induce 
a permanent material modification^^. 

Several theoretical investigations concerning the inter¬ 
action of fs laser pulses with dielectrics can be found in 
the literature. For instance, Brouwer et alJ provided a 
detailed study of fs pulse propagation and electronic ex¬ 
citation during the interaction. Bulgakova et al.^^ evalu¬ 
ated the energy deposition by a single fs pulse, in order to 
model the heat source in subsequent simulations of ther¬ 
mal accumulation. Eaton et al}^ or Jamshidi-Ghaleh et 
al}^, simulated heat accumulation from pulse trains by 
using phenomenological expressions for the deposited en¬ 
ergy obtained from experimental observations. Here, for 
the first time, we confront rigorous quantitative theo¬ 
retical predictions of permanent material modifications 
induced by a train of fs laser pulses with experimental 
results. 

In this paper, we develop a multi-scale model covering 
the key physical processes eventually leading to the ob¬ 
served permanent modification of the glass through heat 
accumulation. This modification may be caused by, e.g., 
structural change, creation of color centers, etc. In our 
model, we simply assume that permanent modification 
occurs for lattice heating above a given threshold tem¬ 
perature. Our approach is based on separation of the 
different time-scales involved, i.e., fs or ps for pulse prop¬ 
agation and energy deposition, and ps or ms for the ther¬ 
mal processes. First, we simulate the highly nonlinear fs 
pulse propagation accounting for multi-photon and col¬ 
lisional ionization as well as the optical Kerr effect. In 
a second step, we extract a map of the deposited en¬ 
ergy density in the focal volume. This energy density 
map is then used as a source term in a heat equation in 
order to compute the ms evolution of the thermal dis¬ 
tribution. The model is evaluated for trains of up to 
500 pulses with microjoule energies and kilohertz RRs 
focused into soda-lime glass to investigate the heat ac¬ 
cumulation regime. Results are directly confronted with 
experiments for quantitative comparison. Good agree- 
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ment is obtained for threshold temperatures located in 
between strain and annealing point of the material. 

Let us start with a more detailed description of fs laser 
pulse propagation in transparent media, thus considering 
the fs-ps timescale. To this end, we resort to an extended 
nonlinear optical Schrodinger equation accounting for 
self-focusing due to the optical Kerr effect and light de- 
focusing via laser-generated conduction electrons^^. For 
the considered case of a linearly polarized light and mod¬ 
erate focusing conditions, the scalar and paraxial approx¬ 
imations are justified. Moreover, by comparing to a more 
complete unidirectional model^^ we checked that higher 
order effects like self-steepening can be neglected. Thus, 
the slowly varying complex optical field envelope ^(r, t, z) 
is governed by the following equation 
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where no and n 2 are the linear and nonlinear refrac¬ 
tive index, respectively, k" corresponds to the second 
order group velocity dispersion in the material, Nc is 
the critical electron density in the conduction band. Eg 
is the band gap, is the laser pulse intensity, 

where the optical field envelope 5 is normalized such 
that I = The time i is the retarded time in a 

frame moving with the group velocity of the pulse at 
center frequency coq and wave number ko = noUJo /c. The 
laser energy absorption by conduction band electrons is 
described by the Drude model with the cross section 
a = koe‘^rc/[nQUJoSome{l E ooqT^)], where Tc is an effec¬ 
tive electron collision time. The evolution of the electron 
density in the conduction band Ne{r, t, z) reads 


diNe = Wpi {Nnt-Ne) + NE-^, (2) 

'^rec 

where N^t is the initial density of valence electrons, Wpi 
is the photo-ionization rate derived from Keldysh for¬ 
mula^, and Tree is a typical recombination time. 

Solving Eqs. (1) and (2) gives the spatial distribution 
of laser energy density absorbed by the CB electrons dur¬ 
ing the whole propagation of a single laser pulse. This 
deposited energy density U{r^z) is assumed to be fully 
transferred to the lattice much faster than the charac¬ 
teristic heat diffusion time. Thus, irradiation by each 
laser pulse increases the lattice temperature T(r, t, z) by 
Tl{t^z) = UjpCv^ where p and Cy are the density and 
heat capacity of the material, respectively. The spatial 
and temporal evolution of T is then determined by 


BT ^ ^ 

— = DV^T+J2TL5{t-jAt). (3) 

j=0 

Here, the first term on the right-hand-side describes the 
temperature diffusion process with the thermal diffusion 


no 

1.5134 

CT [cm^] 

1.05 X 10“^^ 

712 [cm^/W] 

3.86 X 10“^® 

Tc [fs] 

1.83 

Eg [eV] 

3.9 

Tree [fs] 

150 

k" [fs^/cm] 

286 

P [g/cm®] 

2.44 

Nc [cm“^] 

1.05 X 10^^ 

D [cm^/s] 

0.72 

Nnt [cm“^] 

2.1 X 10^^ 

Cv [J/g.K 1 

5.97 X 10"® 


TABLE 1. Material parameters used to simulate femtosecond 
pulse train interaction with soda-lime glass^^“^^. 


coefficient and the last term corresponds to the heat 
source caused by a train of N pulses separated by At. 
Because the expected material modifications are rather 
weak, we neglect any temporal variation of the diffusion 
coefficient D and the source term Tp. Both quantities 
are evaluated at room temperature. 

By solving Eq. (3), the lattice temperature evolution 
of the material irradiated by a fs pulse train is obtained. 
If the thermal diffusion process in Eq. (3) is slower than 
the heating process, the increase in the temperature can 
be significant for long pulse trains containing several hun¬ 
dreds of pulses. We assume a permanent material modi¬ 
fication in regions where the temperature becomes larger 
than a given threshold temperature. This condition de¬ 
fines both the size and the shape of the material modifi¬ 
cation with respect to the laser pulse train parameters. 

Simulations are performed for the soda-lime glass ir¬ 
radiated by trains of 500 pulses with 290 fs duration 
at 1030 nm wavelength and 1.3 pJ single pulse energy. 
Pulse parameters and focusing conditions are chosen cor¬ 
responding to the experiments presented below, where 
the position of the geometrical focus is about 500 pm 
below the sample surface inside the bulk material. All 
relevant material parameters are summarized in Tab. 1. 

Eigure l(a,b) presents spatial and temporal dynamics 
of the single fs pulse propagation in soda-lime glass ob¬ 
tained from Eqs. (1) and (2). The evolution of the beam 
fluence shown in panel (a) confirms a focal spot size of 
a few microns. The on-axis temporal evolution shown in 
panel (b) is typical for focused fs low-energy pulses: a 
short (here ^ 50 fs) intense spike at the pulse front is re¬ 
sponsible for a significant free electron production, which 
in turn defocuses the trailing part of the pulse. 

Erom the single fs pulse propagation simulation, we 
extract the deposited energy density U in order to con¬ 
struct the temperature distribution Tp transferred to the 
lattice, i.e., the source term in Eq. (3). In the simula¬ 
tion, 60 % of the laser energy is absorbed in the ma¬ 
terial, which is close to 55 % measured experimentally. 
The main part (95 %) of the absorbed laser energy is 
deposited in the material located before the focal point, 
even though the absorbed energy density is maximum at 
the focus (see Eig. 1(c)). This particular spatial distribu¬ 
tion is induced by the promotion of electrons to the CB 
during the propagation, because these electrons absorb 
the laser energy very efficiently. After the focal point, 
the intensity rapidly decreases and becomes too low to 
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a: Fluence [J/cm^] b: On-axis intensity [TW/cm^] 
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FIG. 1. (Color online) Fluence (a) and on-axis temporal profile (b) versus propagation distance for a single pulse (1.3 /iJ, 
290 fs at 1030 nm) propagating in soda-lime glass according to Eqs. (1) and (2), for an initial beam diameter of 5 mm and a 
focal length of 20 mm (in air). The deposited energy density by such single pulse is directly proportional to the temperature 
distribution Tl transferred to the lattice (c). The evolved temperature distribution according Eq. (3) after 2 ps, just before 
the second pulse arrives, is shown in (d). The laser pulse propagates from left to right, with geometrical focus at zero. 


create a significant carrier density in the CB. 

Our simulations results of the single fs laser pulse prop¬ 
agation show that the maximum temperature reaches 
T = 360 K at the focus (see Fig. 1(c)), which is much 
smaller than the strain temperature = 800 K or an¬ 
nealing temperature Ta = 830 K in the soda-lime glass^^. 
Therefore, we expect no material modification from a sin¬ 
gle fs laser pulse. However, Fig. 1(d) reveals that 2 ps 
later, i.e., just before the second pulse arrives, the tem¬ 
perature in the interaction zone is still significantly ele¬ 
vated. In fact, the heat transfer in the surrounding ma¬ 
terial is slow, only the spatially very narrow temperature 
peak at the focus has decreased to some extend, and the 
temperature in the overall interaction zone is still well 
above the initial 300 K. It would take more than 100 ps 
to relax to temperatures below 305 K in the whole volume 
after a single pulse interaction. Thus, the interaction vol¬ 
ume will get heated to higher and higher temperatures 
by each pulse of the incident fs pulse train, and a thermal 
accumulation process takes place. 

The lattice temperature obtained after 1 ms (accumu¬ 
lation over 500 pulses at a 500 kHz RR) is shown in Fig. 2. 
The interaction zone is located about 100 pm before the 
focal point. The radial dimension of the heated volume 
in Fig. 2 is much larger than the one affected by the single 
fs pulse (c.f. Fig. 1(c) and (d)). The maximum temper¬ 
ature after 1 ms is close to 1350 K, i.e., well above both 
Tg and Ta^ and a permanent material modification can 
be expected in the region where T > encircled by the 
black line in Fig. 2. 

To compare our numerical predictions with experimen¬ 
tal observations, pulse trains produced by an Yb-doped 
fiber laser (Satsuma, Amplitude Systemes) were focused 
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FIG. 2. (Golor online) Spatial distribution of the temperature 
in soda-lime glass after 1 ms of irradiation by 500 pulses at 
500 kHz RR. The zone with temperature above Ts is encircled 
by the black line. The laser pulses propagate from left to right. 


into a sample of a bulk soda-lime glass with a xlO Mitu- 
toyo objective lens (NA = 0.26; / = 20 mm). Figure 3 
presents transverse images of the interaction volume after 
irradiation with 500 pulses for various RRs. We report no 
measurable modification of the matter for 100 kHz RR; 
permanent modifications of the material clearly appear 
for RRs larger than 200 kHz. The size of the structure 
increases with the RR towards the surface of the sample. 

For a more quantitative comparison, Fig. 4 presents the 
measured and simulated dependence of the size (length 
and diameter) of the material modification versus the 
RR. For a threshold temperature near the strain point 
{Tg = 800 K) or the annealing point {Ta = 830 K), i.e. 
in the glass transition region Tg^ we report good agree¬ 
ment between experiment and model. We did the same 
analysis for different numbers of pulses and pulse ener¬ 
gies, and found the same agreement (not shown here). 
The fact that we find a threshold temperature around 
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model accounts for nonlinear laser pulse propagation as 
well as thermal diffusion. Simulations were performed for 
soda-lime glass, and directly confronted with experimen¬ 
tal results. We report that depending of the RR, thermal 
accumulation in and around the focal volume leads to ma¬ 
terial modifications. It appears that the transition tem¬ 
perature of the glass can be considered as a good thresh¬ 
old for permanent material modifications. The size of the 
simulated modification volume is in good agreement with 
experimental observations in the thermal regime for RRs 
higher than 200 kHz. We expect that our model is appli¬ 
cable to a wide range of dielectrics. Thus, it will provide 
a new baseline for material structuring by fs lasers. 
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FIG. 3. Transverse optical microscopy images of the perma¬ 
nent material modifications induced in soda-lime glass for 500 
pulses for various RRs. 



FIG. 4. Quantitative comparison of experimental (stars) and 
modeling data (squares and triangles) versus RR. The dashed 
lines serve as guides to the eye. 


Tg confirms that the observed modifications are mostly 
due to thermal accumulation, i.e. structural changes and 
related induced stress. For a RR lower than roughly 200 
kHz, the predicted lattice temperature remains below Tg. 
The small modifications observed in the experiments at 
100 and 200 kHz could thus be associated to non-thermal 
processes such as changes in the electronic structure. Ac¬ 
tually, preliminary experimental results for the lowest 
RRs based on Difference Intensity Contrast microscopy 
exhibit slight modifications, corroborating non-thermal 
material modifications in this interaction regime. 

In summary, a multi-scale model has been developed to 
describe permanent material modifications in dielectric 
materials irradiated by femtosecond pulse trains. Our 
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